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"^S ' Abstract 

■^1- . We provide in this article a new proof of the uniqueness of the flow solution to ordinary 

£*\J ' differential equations with BV vector- fields that have divergence in L°° (or in L 1 ), when 

the flow is assumed nearly incompressible (see the text for the definition of this term) . The 
novelty of the proof lies in the fact it does not use the associated transport equation. 
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1 Introduction and statement of our main result 

In 1989, P.-L. Lions and R. DiPerna showed in [TU] the existence and the uniqueness of the 
almost everywhere defined flow solution to an ordinary differential equation of the type: 

VO . for W 1 ' 1 vector fields b with L} oc (M.t,L?f) divergence (along with some technical assumptions). 

tH- ' For such 'singular' vector fields, the only possibility is to solve the equation almost everywhere 

on the space f2 of initial conditions. In that case, one defines a flow X(t,x) : R x il — > O 

CO \ satisfying 

X(t,x) = b(t,X(t,x)) for alii, ,„s 

X(0,x)=x ' l ' 

for almost all £ £ il. An initial time s^O may of course be chosen, and the flow then depends 
parametrically on this initial time s. With a view to simplifying the presentation, we will assume 
henceforth and throughout this article that the field b is time-independent. Our arguments may 
be modified to cover the time-dependent case. 

In the present article, we also adopt a notion of almost everywhere flow solution similar to 
that of DiPerna and Lions. We denote by (X(t, -)#^) (E) — MX{— t,E)) the pushforward of 
the Lebesgue measure A. In the sequel, the vector-field b will always be assumed at least L\ oc . 
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Definition 1 (Almost everywhere, nearly incompressible flows). An almost everywhere flow 
solution to ([2]) is a measurable function X(t,x) : R x SI — > $7 satisfying the following conditions: 

(i) For almost all x £ fl, the map t \-} X(t,x) is a continuous solution to 7 = 6(7) satisfying 
7(0) = a; and 

/or almost all x £ $7, Vi £ R, -X"(t, x) = x + / b(X(s, x)) ds 

Jo 

(u,) -For a?/ i, the measure X(t,-)#\ is absolutely continuous with respect to X, and, for all 
T > 0, there exists a time- dependent function C(T) > such that the following condition 
of near-incompressibility holds 

VT £ R+, V|t| < T, -J—\<X(t,-)#\<C(T)\, 

(Hi) X is a one-parameter transformation group, i.e. satisfies 

X(t, X(s, x)) = X(s + t, x), for almost all x £ 51, Vs, t 

Remark 1. Di Perna and Lions originally define in flOf a flow solution with condition (i) 
replaced by X £ C(R, L 1 ) satisfies the ordinary differential equation ^ in the sense of distribu- 
tion. Their definition is equivalent to ours. It is indeed shown in Jiff] that the original definition 
implies (i), and it can be shown that conditions (ii) and (Hi) together imply that X £ C(R, L 1 ). 

Remark 2. Condition (ii) is in particular satisfied if b is sufficiently smooth, and div(6) £ 
L°° . In that case C(t) = e" dlv ^'"°°' t ' is convenient. Besides this case, for a class of ordinary 
differential equations coming from some particular types of hyperbolic equations, such as the 
Keyfitz-Krantzer system, div(o) is only L 1 , but an estimate of the form (ii) may be established 
using a maximum principle. See the work [6] by Bressan for more details on these systems 
and that by L. Ambrosio, F. Bouchut and C. De Lellis W for a discussion on the relevance of 
condition (ii). 

We now give a brief state-of-the-art survey on the theory of ordinary differential equations 

with vector fields of low regularity. The seminal work |10] by DiPerna and Lions has been 

complemented and extended notably by L. Ambrosio in |2J. Several other authors have made 

important contributions. We would like to specifically cite the work |13| by N. Lerner which has 

inspired our own, present work. To date, the minimal conditions that are known to guarantee 

the existence and the uniqueness of the flow are the BV regularity of the vector field, a L 1 

bound on the divergence together with a near-incompressibility condition (or more classically a 

condition of bounded divergence) of the type (ii). The classical proofs of such results are based 

upon the consideration of the associated transport equation, written either in the conservative 

form 

Fin 

^ + div(6(x) M )=0, (3) 

°-g + b(x) ■ V x u = 0, (4) 



or in the non-conservative form 



both with the initial condition w(0, x) — uq(x). Remark that, for divergence-free fields, the two 
equations coincide. When the existence and the uniqueness of the solution to the transport 
equation is established, for any given initial condition, one deduces the same result for the a.e. 
flow solution to the ordinary differential equation. The key ingredient for the resolution of the 
transport equation is a commutation lemma (first stated in |10j). which states that 

p e * divfbu) — div((6 * Pe) u ) — > in L . 

C. De Lellis and G. Crippa have recently given in [9] a new proof of the existence and uniqueness 
of the flow solution of (J2J , not using the the associated transport equation. Their very interesting 
approach provides regularity estimates for W 1,p vector-fields with p > 1 but seemingly fails for 
W ' vector-fields, unfortunately. Very recently, the approach has been succesfully improved by 
P.-E. Jabin. In [12], this author extends the direct method by C. De Lellis and G. Crippa to the 
case of bounded SBV vector fields (with locally finite jump set for the d — 1 Hausdorff measure) 
in any dimension, and also to two-dimensional BV vector-fields that satisfy a particular, local 
assumption in the direction of the flow. In addition, in [8 , the same author, in collaboration with 
N. Champagnat, has proved that, in the particular case of an ODE corresponding to the Newton 
equation of motion (that is, X = V, V = F(X)), there exists a flow solution if F € H 3 ^ 4 n L°°, 
and this flow is unique among the class of flows obtained by regularization of the problem. 

Before we get to the heart of the matter, some comments are in order regarding our as- 
sumptions, which are slightly different from those of [5]. The important assumption in [5] is the 
bound from above for t > in (ii), in order to prevent contraction of the measure. In addition, 
the group structure (our assumption (iii)) is not required in [5]. In the present state of our 
understanding, such differences seem to us related to the different techniques of proof used in [S] 
and here. The need for (iii) is consistent with our previous work QT] where, in order for the proof 
of uniqueness to be performed, the group structure is necessary in the Lagrangian viewpoint, 
but not in the Eulerian viewpoint. The assumption on the group structure is also present in 
|10j . In the present contribution, we need this assumption because we transform integrals back 
and forth in time (see our argument in Section 2 below). To end this discussion, we also note 
that Assumption (ii) may be slightly weakened if we are only interested by the uniqueness of 
the forward-in-time flow. In that case, we may replace (ii) by 

{ii') ^-_\<X(t,.) # X«X, Vt6[0,T]. 

This is in contrast to the assumptions of [3], which uses a bound from above: X(t, -)#A < C(T) X. 
This owes to the very strategy of our proof, which mostly uses the backward-in-time flow. For 
instance, the most important step ([7]) is performed on the backward flow. In any event, owing 
to the group property, the upper and lower bounds are closely related to one another. The 
two uniform bounds in (ii) are indeed necessary to conclude when the issue considered is the 
uniqueness of both the backward-in-time and the forward-in-time flows. We refer the reader to 
Remark Oat the end of the proof for more details. 

1.1 Main result 

The purpose of this article is to give a new and direct proof of the uniqueness of the a.e. flow 
solution to © for BV vector fields, without arguing on the associated transport equation. We 



adopt the approach already used in [TT] for W 1,1 vector fields. Basically, the commutation lemma 
instrumental in the proof contained in this prior publication is replaced by another strategy of 
proof, namely the introduction of a second variable. This is explained in details in the next 
paragraph. 

Our result is the following: 

Theorem 1. Let b be a BV vector field on the N -dimensional torus T N . If div(b) £ L 1 , then 
there exists at most one a.e. flow solution to @, in the sense of Definition^ 

Remark 3. In fTjj, B. Perthame and I. Capuzzo Dolcetta remarked that the assumption "b £ 
W 1 ' 1 " of the original work by DiPerna and Lions could be replaced by the weaker assumption 
"the symmetric part of Db is a matrix-valued L 1 function". This observation seems to not be 
valid for the present strategy of proof, and more generally in the BV case. The reason is, their 
argument is based on the use of radially symmetric regularization kernels, while the regularization 
kernels we use here for the BV case are typically anisotropic. 

Remark 4. Theorem [7] states uniqueness of the flow. We do not know which condition, besides 
the general assumptions of Theorem [7J makes possible the proof of existence of a flow in the 
sense of Definition [7J using the approach developed here for uniqueness. 

1.2 Main idea of the proof 

To start with, we outline here the proof performed in details in the next section. As already 
said, the proof uses a technique introduced in |11) . In that work, a smooth convolution kernel 
p, with normalized integral, is considered. It is then proved that for any two a.e. flows X and 
Y solutions to the ordinary differential equation with W 1 ' 1 coefficients, 



MfJw>-™&« E T' 



!™-( / / lX(t.x) -Yit.y^pi^O-ldxdp) =0. 



Now, the limit of the integral is 



d 

~dl 



\X{t,x) -Y{t,x)\dx) = 0. 



This shows that, for all t, 



\X(t,x) -Y(t,x)\dx = 0, 



since this quantity vanishes at initial time. The uniqueness of the solution follows. Remark 
that the introduction of the extra- variable y allows to perform the calculation without using the 
transport equation. 

Our aim is to now modify the above approach and treat BV vector fields. For this purpose, we 
use a convolution kernel well adapted to the geometry of the flow and the possible singularities 
of the BV vector field under consideration. In short, we consider the regularization kernel 

—j?p(x, ) with p(x, z) = Fq(\U(x)z\' 2 ) det U(x) , and U(x) = Id + 'yrj( x ) ® w(%) ■ 

e" e 



Here, Fq is a smooth function, 7 is a constant that will be sent to infinity, and r\ is an ap- 
proximation of the direction normal to the jumps of the measure Db. The purpose of such 
a construction is to have a regularization that decreases faster in the direction normal to the 
jumps. The idea of a direction-dependent regularization was first introduced by P.L. Lions in 
|14j . N. Lerner introduced the specific position-dependent regularization used here in |13] with 
a view to simplifying the proof of uniqueness originally given by L. Ambrosio for the BV case. 
His argument, however, is still based upon the equivalence with the transport equation. In the 
present paper, we combine his argument with the approach consisting in introducing a second 
variable, already employed in [TT] for W 1,1 vector fields. 

2 Proof of Theorem [I] 

This section is devoted to the proof of Theorem[TJ We denote by pi{t, ") (resp. ^2 (A •)) tne ^°° 
density of the measure X(—t, -)#A (resp. Y{— t, -)#A) with respect to A. 

Consider now the kernel 

■&!>(*>— )• 

where p is a smooth, compactly supported, function, from T N x T^ to R + which we will make 
precise below. Assume in addition p satisfies J p(x, z) dz = 1 for all x. Our aim is to estimate 

Uff\X {t ,x)-Y { t^ p{ x,^ 



I e( t ) = -Al l\X(t,x)-Y(t,y)\—p(x,—^)p 1 {t,x)p 2 (t,y)dxdy). (5) 



where X and Y are two flow solutions to ©. In the sense of distributions, 

lim/ e (t) = — (f \X(t,x)-Y(t,x)\fn(t,x)fi2(t,x)dx\ . (6) 

This is established using the Lebesgue continuity of the functions Y and /12 at almost every 
point, along with the L°° bound on pi. Remark that the Lebesgue continuity may be used if 
the support of p(x, ■) is not exceedingly stretched in one direction (more specifically, we should 
have some constant c > such that Vx e T N , B(0, c _1 ) C Suppp(x, •) C B(0,c), See [TH] for 
more details). The kernel we shall use satisfies such a condition for all e > 0, even though in 
the limit of a vanishing e, it is infinitely stretched. Our purpose is to show that the limit (J6j) is 

liml e (i) = - / \X(t,x)-Y(t,x)\div(b)(x)p 1 (t,x)p 2 (t,x)dx. 

e^O J 

This will eventually prove the uniqueness of the flow solution to ([2]) using the bounds from 
below on p\ and pi inferred from (ii). To this end, we first perform the change of variable 
(x,y) — > (X(t,x),Y(t,y)) in I e (t), and then differentiate under the integral 

d fff, , 1 , v , ^ X(-t,x)-Y(-t,y). , 

e{) = di( / Is-Plpg^H.s). £ l!J, )dxi 



which we write I e {t) = I^(i) + I^(t), with 

^(*) = - / / l -^d lP (X_ t (x), X - t{x) ~ Y ~ t(y) ) • b(X_ t (x)) dx dy 



,N e 



m = - I I l -^d 2 p(X_ t (x), X - t(x) Y - t{y) )-(b(X. t {x))-b(Y_ t (y)))dxdy 



with the notations X t (x) — X(t, x) and Y t (y) = Y(t, y). Then, we return to the original variables 

(x,y) 

Il(t) = - f f \X t (x)-Y t (y)\^ w d 1 p(x,^-)-b(x)ix l (t,x)n 2 (t,y)dxdy 

Ie(t) = - J f\X t (x)-Y(t,y)\^d 2 p(x,^-)-(b(x)-b(y))n 1 (t,x)n 2 (t,y)dxdy, 

and next use the change of variable z = (y — x)/e 

I e {t) = \X t (x) —Y t (x + ez)\dip(x,z) ■ b(x)fj,i(t,x)ii2(t,x + ez) dxdz (7) 



£(*) = -J J \X t (x) - Y t (x + ez)\ d 2 p(x, z) ■ {b{x) y + £Z)) . . . 

. . . p,i(t, x)p,2(t,x + ez) dxdz. (8) 

We now need to estimate these two terms when e goes to zero. We begin with the easiest of the 
two, namely /^. 

Step 1: Limit of I] 

Because p is smooth, b £ L 1 , and almost all points are Lesbesgue points for the two functions 
Y and /Z2, we can use the Lebesgue dominated convergence theorem and obtain 

toJi W = -/ Wa! )-y t (x)|(/Mx,.)«b).6(*)M 1 (*,x) / i 9 (t >a ;)«to. 

Now 

dip(x, z) dz = — I / p{x 1 z) dz I = , 



since J p(x, z) dz = 1, for all x. Thus, 



lim4 1 (t) = 0. (9) 

e->0 



The treatment for 1^ is more elaborate and will necessitate several steps. 

Step 2: Bound for /f 

We now wish to pass to the limit e — ► in ||SJ. If b were W 1 ' 1 , the limit could easily be 
identified. It would suffice to replace (b(x + ez) — b(x))/e by f Q Db(x + 9ez)- zdO in ((8|), and next 
use the Lebesgue dominated convergence theorem. All this does not require making specific the 
convolution kernel p (See below and [H]). Owing to the presence of the singular part of Db, we 
have to argue more carefully. 

To proceed further, we recall the following result: 

Proposition 1. [from [5, Theorem 1.28, Corollary 1.29]] Let b be a BV vector-field on 

T N . 

(i) The Radon-Nikodym decomposition of its derivative Db writes 

Db = D a b + D s b, with D a b«C N , D s bLC N , 



where the superscript a stand for "absolute continuous part", and s stand for "singular" respec- 
tively. As D a b is absolutely continuous with respect to the Lebesgue measure, we write it 

D a b = d a bdx, 

where d a b is a L 1 matrix-valued fonction. 

(ii) In addition, the polar decomposition of the singular part D s b of the measure Db writes 

D s b = M s \D s b\ , 

where \D s b\ is the total variation of the matrix-valued measure D s b, and M s a matrix-valued 
fonction, such that \M s {x)\ = 1, \D s b\-a.e (the norm used for M s is the norm induced on 
matrices by the Euclidian norm ofM. n ). 

In view of the above decomposition, we now claim that 



lim sup 

£ — >Q 



\X t (x)-Y t (x + ez)\ 



i),,,(.v.z)- ( 6 ( g + g *) Kg) -Ptyx + eOzyz 



pi(t,x)p 2 (t,x + ez) d6 dx dz (10) 



<2C{t) 2 / / \d 2 p(x,z)-M s (x)-z\d\D s b\{x)dz. 
For convenience, we denote by 

ll a = - ((J \X t (x)-Y t (x + ez)\d 2 p(x,z)-d a b(x + e9z)-zd9dxdz, 
in the left-hand side, and 

P s (t)= If \d 2 p{x,z)-M s {x)-z\ d\D s b\(x)dz, 



in the right-hand side. 

To prove our claim, we regularize X, Y, p\ and \x 2 , using some smooth X a , Y a , fj," and 
^2 . Next, we replace (b(x + ez) — b(x))/e by J Q Db(x + 9ez) ■ zd6 (an equality true for almost 
all (x, z)) and perform the change of variable x' = x + eOz (we use it even for the measure Db 
because this is a linear change of variable) . We obtain 



Ie' U (t) 



\X?(x)-Y t a (x + ez)\d 2P (x,z) 



(b(x + ez) -b(x)) 

e 
. . . Hi (t, x)ix 2 (t, x + ez) dx dz 

\X?(x) -Y t a (x + ez)\d 2 p(x,z) ■ Db(x + 6ez) ■ z . . . 

. . . Hi (t, x)n 2 (t, x + ez) dx dz d9 

\X?(x - e9z) - Y t a (x + e(l - 9)z)\ . . . 

. . . d 2 p(x — eOz, z) ■ Db(x) ■ zpi(t, x — e9z)p^(t^x + e(l — 9)z) dxdz d9 . (11) 



Let us decompose if'" in two parts, according to the above Proposition [TJ 

where 

I^(t) = - J J J \X a (t,x-e0z)-Y a (t,x + e(l-O)z)\... 

. . .d2p(x—£0z, z)-d a b(x)-z ^(t,x-e9z)^2(t,x+s(l — 8)z)dxdzd9 
|l£f(t)| < 2C{t) 2 f f f d 2 p{x-e9z,z)-M s (x)-z\D s b\{x)dzd0 

where we have used that \X a — Y a \ < 2 (as we work on the torus). And letting e going to 
zero, we obtain (10) for X a , Y a and the fj,f. Then flTU|) is obtained letting X a ,Y a , nf and fi% 
approximate X, Y, p\ and p 2 , respectively. 

The majoration ([TOj being established, we proceed as follows. Arguing as above for J*, that 
is using the smoothness of p and the fact that almost every point is a Lebesgue point for Y, pi 
and p 2l we obtain 

lim i 2 a (i) = lhx l - I f f \X t (x - e9z) - Y t (x + e(l - 9)z)\ . . . 

. . .d 2 p(x — e$z,z) ■ d a b(x) ■ z pi(t,x — e9)p 2 (t,x + e(l — 0)z) d9 dxdz 
\X t (x) — Y t (x)\d 2 p(x, z) ■ d a b(x) ■ z/J,i(t, x)p 2 (t,x) dxdz 



= \X t (x) -Y t (x)\R a (x)p 1 (t,x)p 2 (t,x)dx, 

with R a {x) — J d 2 p{x, z) ■ d a b{x) ■ zdz. To calculate this term, we integrate by parts and use 
the property \/x € T , J p(x, z) dz = 1 

= -Y^{x) f p(x lZ )dz = -div a b. (12) 

So we have obtained 

limi £ 2 a = [\X t (x)-Y t (x)\div a b(x)p 1 (t,x)p 2 (t,x)dx. (13) 






The next step consists in proving that the right-hand side of (|10j) may be chosen arbitrarily 
small. 

Step 3: A bound on the singular part 

In order to estimate the right hand side of (|10[) . we now use a geometric information, namely, 
the special form of M s (x), proved by G. Alberti pQ. 



Theorem 2. [Alberti's rank one Theorem, [5j Theorem 3.94]] Let b be a BV vector-field 
defined on T N , and write Db = D s b + D a b the Radon- Nikodym decomposition of its gradient. 
Consider D s b = M s \D s b\ the polar decomposition of the singular part as in Proposition]]^ Then, 
M s is of rank one \D s b\-almost everywhere, that is, there exists two vector-valued functions £(, 
and r\b, both \D s b\-measurables, such that £b and r\b are unit vectors \D s b\-a.e. and satisfy 

M s (x) = £b{%) ® Vbix), \D s b\ — almost everywhere, 

where ^ ® r\b denotes the linear map z t-¥ {r)b, z)£i>- 

Corollary 1. As a consequence, the singular part of the divergence is 

div s 6= (£,r))\D s b\. 

If we assume that the divergence of b belongs to L , it follows that 

(£,T)) = 0, \D s b\ — almost everywhere, 

a property that will be crucial in the sequel. 

Using the decomposition provided by Theorem 2, we rewrite our bound in (|10p . which we 
denote I 2 (t) in the sequel 

P s {t)<2C{tf J J\{d 2P {x,z),ib{x))\\{r 1 b{x),z)\d\D s b\{x)dz. 

In order to render the right-hand side arbitrarily small, we now make specific our convolution 
kernel p. We choose 

p{x,z)=F a (\U(x)z\ 2 )det(U{x)), 

where Fo is a smooth, compactly supported, non negative function such that 
J R iv Fo(\z\ 2 )dz — 1, and U is a smooth, matrix- valued function, such that U(x) is an orientation 
preserving matrix for all x. Note that owing to the presence of the determinant, the integral 
of p(x, •) remains equals to one independently of x. The dilation matrix U{x) is set to U(x) = 
Id + 777(2;) CS> f){x) (with the notation a (g> b for the endomorphism x — ► (b,x)a), where r\ is a 
smooth vector-valued function. On the jump of the measure Db, r\ will be chosen later as an 
approximation of the direction normal to the jump set. The factor 7 will be chosen as large as 
possible. It may possibly depend upon x and be large only on a neighbourhood of the singular 
set of the measure Db, but we for simplicity of the calculation we will not use that not essential 
possibility here. 

The partial derivative of p writes 

d 2 p{x,z) = 2F^(\U(x)z\ 2 ){U(x)z,U(x)-)det{U(x)). 

We use this in the bound on I 2 (t) to obtain 

P s it) <C f f\F( ) (\U(x)z\ 2 )\\(U(x)z,U(x)^(x))\(f ]b (x),z)\det(U(x))d\D s b\(x)dz, 



where here and below C denotes various irrelevant constants. To simplify this term, we perform 
the change of variable z — > U(x)z, and obtain 

i 2 M<C I [\F^\z\ 2 )\\(z,U{x)Um\{v b {x),U- l z)\d\D s b\{x)dz. (14) 



We next intend to use the special form U(x) = Id+'y(x)r)(x) ®r)(x) , to bound from above the 
two scalar products. Let us first formally illustrate our argument, performing our calculation 
with r\ = rjb, as if Vb were smooth. In this case, 

\(z,U(x)£ b (x))\ = \(z,t; b (x)}\<\z\, 

because {vb,£,b} = and £b has unit norm. For the second scalar product, 



\{ Vb (x),U- 1 (x)z}\ 



1 + 7 



(r) b (x),z)\ < 



1 + 7' 



because U 1 = Id 



7 



1 + 7 
m 

i 2 s (t) < 



Vb <x> Vb- 



Inserting these bounds in (ITi|) . we obtain 

C 



1 + 7 



\F'(\zn\D s b\(x)dz< 



C(F ,b) 

1 + 7 : 



where the constant C(Fo,b) depends only of Fq and 6. It remains then to let 7 to infinity to 
obtain I^{t) — and conclude our (formal) proof. 

We now modify the above formal argument using an approximation 77 of rjb, instead of rjb 
itself. First, we remark 

\{z,U{x)tib)\ = |(*,6+7<6,»7to>l 

< (1 + 7KMI) \z\ 

< 0- + -y\(£b,v-Vb)\)\z\ 

< (i + 7 h-%l)kl, (is) 

where we have used (£b,Vb) — 0, |-D s 6|-a.e (from Corollary [1]), and that £&, rjb are unit vectors. 
To bound the scalar product \(r)b(x), U z)\, we decompose z in z = z„ + z±, where z„ is the 
projection of z on M rj 



\{Vb,U- l z)\ 



Vb,z 



1 + 7 



(V,z}v 



Vb, Z± + 



1 i V 

1 + 7 



< 



< 



Vb ~V,Zx + 



\Vb-v\ + 



1 + 7 



V,z± + 



1 + 7 



1 + 7 



(16) 



From (fT5|) and ([T^j) we deduce 

|^,C/(^6>||<^,C/-^)|< ( 2|77 — r? & | 



1+7 



+ l\v~Vb\ 2 ) \z? 



10 



We insert this bound in (IT?)) and obtain 

P s {t) < Cj ^\ v -r ]b \ + T ^+ 1 \r 1 - Vb \^jFi(\z\ 2 )\zfdz)j\D s b\( X ) 

~ C, ( i? o)( T ^ + (l + 2 7 )|h-r /b ||^6|( a; )), (17) 

because the integral J Fq(\z\ 2 )\z\ 2 dz is fixed, and both r\ b and r\ are unit vectors. 
We finally show that 

inf C T i_ + (l + 2 7 ) [\V-Vb\ \D s b\(x))=0. (18) 

7>0,?7 smooth ^ 1 + 7 J I 

To this end, we first choose 7 such that 1/(1 + 7) i s small, and then construct a smooth function 
r], sufficiently close to rjb on the support of D s b so that (1 + 27) J \rj — r]b\ \D s b\(x) is also small 
(use for that classical approximation theorem with respect to the Radon measure D s b). Note 
that r\ can be arbitrarily extended to the whole torus as its value outside the support of D s b is 
irrelevant. This concludes the proof of the convergence of the right-hand side of (fT0"|) to zero. 

Step 4: Conclusion Collecting all the previous results, we obtain 

— I \X t {x) — Y t {x)\n\(t,x)n2(t,x)dx = 

(19) 
\X t (x) - It (a:) I div(b)(x)ni(t, x)[i 2 (t, x) dx. 

where we have replaced div a (6) by div(&), since we are dealing with vector fields b having at 
least, divergence in L 1 . If div(fe) G L°°, then 

— / \X t (x) -Y t {x)\n 1 (t,x)n 2 (t,x)dx < C I \X t (x)~Y t (x)\fj, 1 (t,x)fj, 2 (t,x)dx. 

Since the integral in the right hand side vanishes initially, we conclude that 

\X t (x) - Y t (x)\iii(t, x)n2(t, x) dx = 



and finally that X(t, •) = Y(t, ■) a.e. in x since the (Xi are bounded away from 0. Note that, as 
usual, if only the solution at positive times if of interest, an assumption on the negative part 
div(6)~ of the divergence suffices to conclude. 

When only the weaker hypothesis div(6) g L 1 holds, we have to slightly adapt the above 
argument. We choose a smooth compactly supported function 4>{x), insert a factor (j)(X(t, x)) 
in the integral ((5j) defining I e . We now estimate 

#(*) = jAJ J <t>{Xt{x))\X t {x) -Y t (y)\^ w p(x,^-)n 1 (t,x)n 2 (t,y)dxdy^ . 



The above argument carries over to the present case. An equality similar to (1191) is obtained 

— U(X t (x))\X t (x) -Y t (x)\m(t,x)n2(t,x)dx 

' (20) 

4>(X t (x))\X t (x)-Y t (x)\div(b)(x)iJ,i(t,x)n2(t,x)dx 
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which can also be written (using the change of variable x — X t (x)) 

^ [<t>{x)\x-Y t {X_ t {x))\ii2{t,X_ t {x))dx 

J (21) 

(x)\x-Y t (X^ f (x))\dw(b)(X- t (x))Li2(t,X- t (x))dx. 

We next dehne u(t,x) = \x — Yt(X^t(x))\n2{t,X-t(x)). Equation (|2Tj) holding for all <fi, it 
follows that 

— +dw(b)(X_ t (x))u = 0, (22) 

at 

in the distributional sense. There is no derivative of u with respect to x in the equation, so that 
the variable x is only a parameter. Since div(6) 6 L 1 and condition (ii) holds, we have 

T 

\ div (b)(X(-t,x))\dtdx < +oo 
for all time T. So that, for almost all x, L \ div (b)(X(—t,x))\dt < +oo. Therefore equation 



(|22|) is well-posed for almost all x, and since by construction its solution u vanishes at initial 
time, it vanishes for all time: u(t,x) = for all t, a.e. in x. This concludes the proof: X = Y. 

Remark 5. Given the above argument, the reader may now understand that, when the unique- 
ness of only the flow for positive times is under study, the upper bound X(t, -)#A < C(T)X, 
for times t > in (ii) of Definition [7] may be somehow relaxed. In that case, and as briefly 
announced in the introduction, we may only assume that 

(n') ^yA<X(i,-) # A«A, Vte[0,T]. 

where the symbol << means here absolutely continuous with respect to. The group property 
allows to equivalently state (composing the previous inequality with X(—t,-)) that 

(ii') A « X(-t, -) # A < C(T)A, Vt e [0, T] . 

In that case, we have 

< fJ>2(t, x) < C(T) , a.e. in x, 

and also < /i2(t, X(—t,x)) almost everywhere in x since, considering 

(ii'), \(X(t,A))<C(T)X(A) 

for all mesurable set A. This suffices to show that, for all t > 0, and almost everywhere in x, 
u(t,x) = and thus conclude X(t,x) — Y(t,x) in the above argument. 

Remark 6. As pointed out by L. Ambrosio in J2j/, the use of the Alberti rank one theorem can 
be circumvented. In our proof, it is possible to use the following, much simpler ingredient. For 
any matrix M , we have 

inf I ' \(Mz,Wp(z))\dz= |Tr(M)|, 
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where the infimum is taken over all smooth kernels with total mass one. Essentially applying this 
result to the matrices D s b{x) we manipulate in the proof, and using Tr(D s b(x)) = drv(b")(x) = 0, 
we may obtain an estimate analogous to (|17[) . thus (|18p. We then conclude our argument 
similarly. We however believe that considering the Alberti rank one Theorem helps to better 
understand the geometry of the problem. 

Acknowledgements. 

This work was mostly completed while the second author was visiting the Institute for Mathe- 
matics and its Applications and the Department of Mathematics of the University of Minnesota. 
The hospitality of these institutions is gratefully acknowledged. Both authors would also like to 
thank Pierre-Emmanuel Jabin for helpful discussions. 

References 

Giovanni Alberti. Rank one property for derivatives of functions with bounded variation. 
Proc. Roy. Soc. Edinburgh Sect. A, 123(2) :239-274, 1993. 

Luigi Ambrosio. Transport equation and Cauchy problem for BV vector fields. Invent. 
Math., 158(2):227-260, 2004. 



[3 



[6 



[9 

[io; 



Luigi Ambrosio. Transport equation and Cauchy problem for non-smooth vector fields. 
Lecture Notes of the CIME Summer school in Cetrary, June 27- July 2, 2005, available at 
http://cvgmt.sns.it/papers/amb05/Cetrarotext.pdf, 2005. 

Luigi Ambrosio, Francois Bouchut, and Camillo De Lellis. Well-posedness for a class of 
hyperbolic systems of conservation laws in several space dimensions. Comm. Partial Dif- 
ferential Equations, 29(9-10):1635-1651, 2004. 

Luigi Ambrosio, Nicola Fusco, and Diego Pallara. Functions of bounded variation and free 
discontinuity problems. Oxford Mathematical Monographs. The Clarendon Press Oxford 
University Press, New York, 2000. 

Alberto Bressan. An ill posed Cauchy problem for a hyperbolic system in two space di- 
mensions. Rend. Sem. Mat. Univ. Padova, 110:103-117, 2003. 

Italo Capuzzo Dolcetta and Benoit Perthame. On some analogy between different ap- 
proaches to first order PDE's with nonsmooth coefficients. Adv. Math. Sci. Appi, 6(2):689- 
703, 1996. 

Nicolas Champagnat and Pierre-Emmanuel Jabin. Well-posedness in any dimension for 
hamiltonian flows with non BV force terms. Accepted for publication in Comm. P.D.E., 
available at http://arxiv.org/abs/0904.1119, 2009. 

Gianluca Crippa and Camillo De Lellis. Estimates and regularity results for the DiPerna- 
Lions flow. J. Reine Angew. Math., 616:15-46, 2008. 

Ronald J. DiPerna and Pierre-Louis Lions. Ordinary differential equations, transport theory 
and Sobolev spaces. Invent. Math., 98(3):511-547, 1989. 



13 



[11] Maxinic Hauray, Claude Le Bris, and Pierre-Louis Lions. Deux rcmarques sur les nots 
generalises d'equations differentielles ordinaires. C. R. Acad. Sci. Paris, 344:759-764, 2007. 

[12] Pierre-Emmanuel Jabin. Differential equations with singular fields. Available at 
http: //matril .unice . f r/^jabin/odef inal . pdf , 2009. 

[13] Nicolas Lerner. Transport equations with partially BV velocities. Ann. Sc. Norm. Super. 
Pisa CI. Sci. (5), 3(4):681-703, 2004. 

[14] Pierre-Louis Lions. Sur les equations differentielles ordinaires et les equations dc transport. 
C. R. Acad. Sci. Paris Ser. I Math., 326(7):833-838, 1998. 

[15] Elias M. Stein. Singular integrals and differentiability properties of functions. Princeton 
Mathematical Series, No. 30. Princeton University Press, Princeton, N.J., 1970. 



14 



